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Fitting a distribution to Value-at-Risk and 
Expected Shortfall, with an application to 

covered bonds 

Dirk Tasche* 


Covered bonds are a specific example of senior secured debt. If the issuer of the bonds 
defaults the proceeds of the assets in the cover pool are used for their debt service. If in 
this situation the cover pool proceeds do not suffice for the debt service, the creditors of the 
bonds have recourse to the issuer’s assets and their claims are pari passu with the claims of 
the creditors of senior unsecured debt. Historically, covered bonds have been very safe invest¬ 
ments. During their more than two hundred years of existence, investors never suffered losses 
due to missed payments from covered bonds. From a risk management perspective, therefore 
modelling covered bonds losses is mainly of interest for estimating the impact that the asset 
encumbrance by the cover pool has on the loss characteristics of the issuer’s senior unsecured 
debt. We explore one-period structural modelling approaches for covered bonds and senior un¬ 
secured debt losses with one and two asset value variables respectively. Obviously, two-assets 
models with separate values of the cover pool and the issuer’s remaining portfolio allow for 
more realistic modelling. However, we demonstrate that exact calibration of such models may 
be impossible. We also investigate a one-asset model in which the riskiness of the cover pool 
is reflected by a risk-based adjustment of the encumbrance ratio of the issuer’s assets. 

Keywords: Covered bonds, expected loss, asset value distribution, quantile, expected short¬ 
fall, method of moments, two-parameter distribution family. 


1. Introduction 

Covered bonds are special debt instruments 1 that are normally served by their issuer. But if the issuer 
defaults and is no longer able to make the payments due to the covered bonds holders, the debt service 
is continued on the basis of the proceeds of the assets in the pool covering the bonds. If in the case of 
the bond issuer’s default the proceeds of the cover pool assets do not suffice for the debt service and the 
bond holders claims are not met by the liquidation of the assets, the bond holders still can take recourse 
to the remaining assets of the issuer. Their claims are then pari passu with the claims of the holders of 
the issuer’s senior unsecured debt. 
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Thus the covered bonds holders have a privileged position compared to most other of the issuer’s creditors. 
Indeed, no loss due to missed payments was ever observed for European covered bonds (EPRS, 2015). 
While this is good news for covered bonds investors, it makes estimation of expected loss for covered 
bonds difficult. 

Covered bonds have been popular investments for a long time especially in Europe. As a consequence of 
the recent financial crisis, covered bonds became an even more important funding instrument for banks, as 
such replacing instruments like securitised products. Therefore, asset encumbrance in the banks’ balance 
sheets due to the higher demand of collateral is increasing. This has led to concerns that expected loss for 
creditors of senior unsecured debt is also increasing (CGFS, 2013). But it is not clear how to efficiently 
quantify this effect for the purpose of loss provisioning or stress testing. 

This paper draws on recent work on fitting lognormal asset value distributions to given probability of 
default (PD) and loss-given-default (LCD) estimates (Yang, 2015) and on asset encumbrance (Chan- 
Lau and Oura, 2014) to suggest simple and intuitive models for covered bonds that allows quantitative 
assessment of expected loss and the impact of asset encumbrance. 

We find that one-period lognormal two-assets model in some situations cannot exactly be calibrated even 
if different time horizons for bond issuer and covered bond holders are taken into account. The proposed 
adjusted lognormal one-asset model represents an acceptable workaround for the calibration issues, but 
only as long as the credit quality fo the cover pool is only moderately better than the bond issuer’s credit 
quality. 

The paper is organised as follows: 

• In Section 2, we revisit the PD and LGD fitting observation by Yang (2015), proving that there is 
a unique solution for the two model parameters and placing the approach in the general context of 
the method of moments. 

• In Section 3, we extend the lognormal one-asset model by Chan-Lau and Oura (2014) to a lognormal 
two-assets model and provide formulae for efficient numerical computation of expected losses for 
covered bonds, senior unsecured debt and junior debt with this model. 

• Computations for the two-assets model are straightforward when the model parameters are known. 
However, we show in Section 4 that the calibration of the parameters leads to a non-linear equation 
system that may have no exact solution at all. 

• In Section 5, we revisit the one-asset model by Chan-Lau and Oura (2014). We adjust the model 
for portfolio heterogeneity in terms of risk by modifying the asset encumbrance ratio related to the 
covered bonds to reflect the expected loss of the cover pool. We show that in the case of a homoge¬ 
neous portfolio the adjusted lognormal one-asset model is equivalent to a comonotonic special case 
of the two-assets model discussed in Section 3. Unsurprisingly, therefore, exact calibration of the 
model again is not always possible. 

• In Section 6, we provide some numerical examples illustrating the impact of asset dependence, of 
asset encumbrance, and of approximating the two-assets model with the adjusted one-asset model. 

• The paper concludes with some comments in Section 7. 

• Appendix A provides background information on mean-variance matching. 
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2. Quantile-Expected Shortfall matching 


The Method of Moments 2 probably is the most popular approach to fitting a distribution to a finite set 
of given characteristics. But it is not the only conceivable approach. Hosking (1992), for instance, argued 
that matching with so-called L-moments might provide results that are better aligned with the results 
of standard goodness-of-fit tests. In a finance-related context, Tasche (2009) suggested quasi-moment 
matching for the probability-of-default curves of rating or scoring models in order to appropriately reflect 
the discriminatory power of the rating models. In a recent paper, Yang (2015) observed that lognormal 
distributions may be uniquely determined by specified PD (probability of default) and recovery rate (or, 
equivalently, LGD). 

Assume that A is an integrable real-valued random variable which is interpreted as showing the value 
of a financial asset at some future point in time - say one year in the future. Of course, the asset value 
interpretation holds only if A is a non-negative random variable. Nonetheless, the following concepts 
apply to any real random variable X: 

• Let D £ K. (interpreted as debt due for repayment at the time horizon). Then the PD (probability 
of default) associated with D is 

PD = P[A < D). (2.1a) 


• Let D e R\{0} such that PD > 0. Then the RR (recovery rate) associated with D is 3 


RR 


E[X | X <D\ 
D 


1 {.Y<U}] 

DP [A < D] ' 


(2.1b) 


Stating the recovery rate is equivalent to stating the LGD (loss given default) which is defined as LGD = 
100% - RR. 


Remark 2.1 In case of X having a distribution with positive density, by (2.1a) the threshold value D is 
the quantile at level PD of X, i.e. 

D = qpo{X), 

with q a defined in general by q a (Z ) = min{x : P [Z < z] > a} for 0 < a < 1 and any real-valued random 
variable Z. In the finance community, quantiles associated with loss variables are often referred to as 
Value-at-Risk (VaR). 

Similarly, in case of X with a distribution with positive density, the recovery rate associated with D is 
closely related to the expected shortfall (ES) of X at level PD: 

RR = ESpg(A) ^ 

with ES q defined in general by 

not 

ES a (Z) = a -1 / q u {Z) du 
Jo 

= E[ZjZ< q a (Z )] - (q a (Z) - E[Z \ Z < q a (Z )]) { ^ [Z ~J* a(Z)] - l) . 

for 0 < a < 1 and any real-valued random variable Z. IfP[Z < q a (Z)] = a, it follows that 

ES a (Z) = E[Z\Z <q a {Z)\, 

2 See Appendix A for a brief description of its simplest version, the mean-variance matching. 

J The indicator function Is of the event E is defined by 1 /.;(T) = 1 for u £ £ and 1 (Y) = 0 for u (j E. 
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while in general we only have ES a (Z) < E[Z | Z < g a (Z)]. In the finance community, ES often is 
associated with loss variables and therefore defined as conditional expectation with respect to the upper 
(right-hand side) tail of the distribution (Acerbi and Tasche, 2002). In the literature, ES is referred to 
also as Conditional VaR (Rockafellar and Uryasev, 2002) or Average VaR (Follmer and Schied, 2011). 

In the following, we will study the problem of how to fit a distribution to given PD and RR (or equivalently 
to given quantile and ES). Hence, we will have two equations to determine the distribution. That is why 
is is natural to consider two parameter location-scale distribution families which are defined as follows. 

Definition 2.2 Let X be a real-valued random variable with given distribution. Then the set of the 
distributions of the random variables m + sX, m £ R, s > 0 is called the location-scale distribution 
family associated with the generating variable X. 

By Remark 2.1, all observations we make in the following for quantile-ES matching of two-parameter 
distribution families immediately apply also to PD-LGD matching. 


2.1. Quantile-ES matching with location-scale distribution families 


Assume that values 0 < a < 1 for the quantile confidence level, q £ R for the quantile, t < q for 
the ES as well as a location-scale distribution family with generating random variable X are given. In 
addition, we assume that X is integrable such that ES a (X) is well-defined and finite. We also assume 
that ES a (X) < q a {X) holds. 

If Y = m + s X with mel and s > 0 is a representation of any element of the location-scale family then 
its a-quantile and ES Q respectively are given by 


Qa(Y) = m + sq a (X), 
ES a (Y) = m + sES a {X). 

Hence, solving q = q a (Y) and t = ES a (Y ) for m and s gives 

q - t. 


m = q — 


q a (X) - ES a {XY 
q-t 


q a (X). 


q a (X) - ES a (X) 

In the special case where X is standard normal with quantile 4 $ _1 (a) = q a (X), we have 5 

(a)) 


ES a {X) = -- 


< 4 (a). 


On the one hand, this provides a proof of the well-known inequality 

<h(— a) < , for a > 0. 

a 

On the other hand, it follows that in the standard normal case (2.2b) reads as follows: 

a(q-t) 


s = 


a<h 1 (a) + <^( < I > 1 (a))’ 
at <h -1 (o:) + q (/?(<l>~ 1 (a)) 
a < t)~ 1 (a) + <^($ _1 (a)) 


4| t> denotes the standard normal distribution function, 
denotes the standard normal density. 


(2.2a) 


(2.2b) 


(2.3) 


(2.4) 
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Since the location parameter of a location-scale distribution family can be any real number, negative 
numbers may be included in the support sets of some or all distributions in the family. Therefore, in 
general location-scale distribution families cannot be used for fitting distributions on the positive real 
axis or the unit interval. 


2.2. Quantile-ES matching with lognormal distributions 


A simple way to generate two-parameter families of distributions on the positive real half-axis is to take 
the exponential of a location-scale family in the sense of Definition 2.2. Hence we study distributions of 
the shape 

Y = exp(m + sA), (2.5) 

where A denotes the random variable generating the location-scale family and we assume in £ B and 
s > 0. 


We want to fit a lognormally distributed random variable Y, as given by (2.5) with X standard normal, 
to given values of 0 < a < 1 for the confidence level, q for the a-quantile and t < q for the ES a . Since 
Y is positive we must also require 0 < t < q. We then have to solve the following non-linear equation 
system for to and s: 

q = q a (exp(m + s A)) = exp(m + s <f> -1 (a)), 

, , . (2.6a) 

t = ES a (exp(m + s A)) = aT 1 exp(m + s 2 /2) ° s ^ m — s). 

See, e.g., Yang (2015, Appendix A) for a derivation of the second row of (2.6a). Rearranging the first row 
of (2.6a) for m and substituting the resulting term for to in the second row gives 


m = log(q) — s<& 1 (a), 

0 = g$($ _1 (a) — s) — at exp(s <J> _1 (o!) — s 2 /2). 


(2.6b) 


Proposition 2.3 For a, b £ R. fixed define the function f : R —> R by 

f(s) = $(a — s) — b exp(a s — s 2 /2). 

If 0 < b < ^(a) then there is exactly one Sq £ R such that f(so ) = 
inequality 


0 < s 0 < 


V{a) 

b 


+ a. 


0 . 


This number sq satisfies the 


Proof. Observe that / has the following properties: 

(i) /(0) = $(o) - b > 0. 

(ii) lim^oo /(s) = 0. 

(iii) / is continuously differentiable with 

f(s) = exp(as — s 2 /2) (bs — ab — ip(a)). 


By (iii) the function / has a unique global minimum at si = + a. By (ii) it must then hold that 

f(s i) < 0. By continuity of / and (i), there is So € (0, si) such that /(so) = 0. This zero of / is unique 
because / has got only one extremum. □ 


By Proposition 2.3, for each triple (a,q,t) with 0 < a < 1 and 0 < t < q there is exactly one solution 
(to, s) € K x (0, oo) of (2.6b) and therefore also (2.6a). The component s of this solution satisfies the 
inequality 


0 < s < 


1 ( a ))g 

at 


$ : (a). 


(2.7) 
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The left-hand side of inequality (2.7) can be somewhat sharpened by making use of the fact that by the 
definition of ES a and the first equation of (2.6b) we have 

t = ES a (exp(m + sX)) < E[exp(m + s A')] = e m+s2/2 = qe 8 " / 2 -* _1 (“) s . 

Some algebra shows that in the case of ($~ 1 (a)) 2 > 2 log (q/t) this implies 

s > $ _1 (a) + \J ($ -1 (a)) 2 - 2 log (q/t). (2.8) 

In the case of ($ -1 (a)) 2 < 2 log(g/t) there is no further constraint for s, i.e. we just have s > 0. 

Another sharper lower bound for s (under a condition on the given quantile q and ES t) can be concluded 
from the following theorem which formalises an observation made by Das and Stein (2013) in the context 
of securitisation. 


Theorem 2.4 Let Z be real random variable and q > 0 such that P[Z < 0] = 0, P [Z < q] > 0 and 
z K>• P[Z < *] is convex for z £ (0, q]. Then it holds that 


nz\z<q] > |. 


Proof. Observe that 

E[Z 1 {Z<q}] 

For any 0 < r < q and r < z < q we have 

P [Z < z] < 


= qP[Z < q] — / P [Z < z\ dz. 


Jo 

by the convexity of * i—>• P[Z < z] that 

«^ P[ Z<r] + ^P[Z<q]. 
q—r q—r 


By r —>• 0 this implies for 0 < z < q that 


(2.9) 


P [Z <z] < ^ P [Z < q}. 

Making use of this inequality in (2.9) gives 

E[Zl {z < q} \ > qP[Z<q]- qP[Z 2 - q] 

_ qP[Z<q] 

2 

Dividing by P [Z < q] now proves the assertion. □ 

Assume that in (2.6a) we have t < q/2. Because lognormal densities are uni-modal Theorem 2.4 then 
implies that the mode of the distribution is smaller than q , i.e. 


R m -s 2 


< q- 


Making use of the first equation of (2.6b) this gives 

qe -*-\a.)s- s 2 K q ^ a > -$-!(a). (2.10) 

For a < 0.001 hence s will necessarily be greater than 3 if the condition t < q/2 on the input data holds. 

More generally, in terms of recovery rates and LGDs, Theorem 2.4 suggests that for ’intuitive’ asset 
value distributions with unimodal densities the expected recovery rate RR is not smaller than 50% (or, 
equivalently, the LGD is not greater than 50%) as long as the default threshold D is not greater than the 
mode of the asset value density. 
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3. A simple model of covered bonds losses 


There is not much literature on mathematical pricing or credit loss models for covered bonds. 

• A detailed pricing model based on time-continuous modelling of defaults and replacements of de¬ 
faulted assets in the cover pool is presented by Kenyon (2009). 

• A one-period model for recovery rates from a bank portfolio consisting of one pool of encumbered 
assets and one pool of unencumbered assets is discussed by Chan-Lau and Oura (2014, Section II D). 
Chan-Lau and Oura model the value of the whole bank portfolio by one random variable. This 
implies an assumption of homogeneity between the encumbered and the unencumbered pools such 
that both sub-portfolios can be modelled with a single asset value variable. 

Driven by the problems with calibration of the model in the absence of proper calibration data (i.e. 
default and loss observations for covered bonds), we follow here the simple approach of Chan-Lau and 
Oura (2014) but introduce several refinements: 

• We start in this section with a model with two asset value variables: one for the encumbered pool 
(e.g. the cover pool for covered bonds) and one for the unencumbered pool, i.e. the part of the 
portfolio which complements the cover pool. 

• We then assume that the bank’s total debt consists of three components: secured (by the cover pool, 
overcollateralised) debt, senior unsecured debt, and junior unsecured debt. Covered bonds usually 
are senior secured debt: In case of the cover pool being insufficient for the service of the covered 
bonds, the uncovered claim from the covered bonds debt is pari passu with the senior unsecured 
debt and higher-ranking than the junior unsecured debt. 

• In addition, in Section 5 we refine the one-asset model of Chan-Lau and Oura (2014) by an adjust¬ 
ment of the encumbrance ratio in order to take into account the risk profile of the cover pool and 
also study the three components debt structure described above in the two-assets case. 

3.1. One-period, two assets approach 

In this paper, we only study the case where both senior and junior (i.e. subordinated) debt are present 
in the portfolio of the covered bonds issuer. From the formulae of that case, the formulae for the case 
without junior debt are readily obtained by setting S = U and U = 0 in the following equations. In 
mathematical terms, the setting is as follows: 

• C is the amount of senior secured (by the cover pool) debt in the issuer’s balance sheet. S is the 
amount of senior unsecured debt. U is the amount of remaining debt, assumed to be subordinated 
unsecured debt. 

• v is the level of over-collateralization of the covered bonds, v = 20% means that the value of the 
collateral is 120% of the face value of the covered bonds. For regulated covered bonds, in general 
there is a minimum requirement of over-collateralisation. Below, the parameter v is used for model 
calibration. 

• X is the future value of the encumbered cover pool assets (securing the covered bonds), Y is the 
future value of the remaining assets, and Z = X + Y is the future total value of the issuer’s portfolio. 
As asset values, X and Y cannot be negative. We consider the values of X , Y and Z as unpredictable 
and, therefore, treat them as random variables. 

• Rc is the recovery rate of the covered bonds, Rs is the recovery rate of the senior unsecured debt, 
Ru is the recovery rate of the junior unsecured debt. Rc, Rs and Rjj are random variables. 

• Lc = 1 — Rc is the loss rate of the covered bonds, Ls = 1 — Rs is the loss rate of the senior 
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unsecured debt, Ljj = 1 — Rjj is the loss rate of the junior unsecured debt. 

In the simplifying one-period context, the following events at the end of the observation period (described 
in mathematical terms) cause financial losses to the creditors of the junior debt, the senior unsecured 
debt or the covered bonds: 

• The issuer defaults, but the total assets suffice to pay for all senior debt (covered bonds and 
unsecured senior debt): 

z <c + s + u, z>c + s. 

The loss rates experienced by the different classes of creditors then are as follows: 

La = 1 - Z ~ ( 0 + S) , La = L s = 0. (3.1a) 


• The issuer defaults, the total assets do not suffice to pay for all senior debt (covered bonds and 
unsecured senior debt). The cover pool suffices to pay for the debt service of covered bonds. We 
assume that the surplus from the cover pool (i.e. the amount of cover pool assets not needed to pay 
for the covered bonds) is given to the creditors of the senior unsecured debt (which, in practice, 
might only happen after many years). This scenario can be described as 

Z <C + S + U, Z <C + S, X > C. 

The implied loss rate are as follows: 

Lc = 0, Ls = 1-——, Lit = 1. (3.1b) 


• The issuer defaults, the total assets do not suffice to pay for all senior debt (covered bonds and 
unsecured senior debt). The cover pool does not suffice to pay for the debt service of the covered 
bonds: 

Z<C + S + U , Z<C + S, X < C. (3.1c) 

In this scenario the creditors experience these loss rates (Chan-Lau and Oura, 2014, see Figure 1): 


Lc 


{C -X){S + C- Z) 
(S + C-X)C 


Ls 


S + C-Z 
S + C-X ’ u 


1 . 


(3-ld) 


At first glance, the case where in total the value of the issuer’s assets is sufficient to pay all debt, but due 
to the asset encumbrance by the covered bonds there is not enough liquidity to pay for the unsecured 
debt (Z > C + S + U, X > C,Y < S + U) also might be considered a loss event. For the purpose of this 
paper, however, we do not treat this event as a loss event because in the longer term the creditors of the 
unsecured debt should be fully paid from the surplus of the cover pool. 

Similarly, the event Z > C+S+U , X < C is not a loss event either because the covered bonds holders’ first 
claim is against the bond issuer which is solvent although the cover pool (despite over-collateralisation) 
would be insufficient to pay for the covered bonds. 


Remark 3.1 An important consequence of the loss assumptions (3.1a), (3.1b), and (3.Id) is that the 
bond issuer suffers a loss (i.e. the issuer defaults) if and only if 

Z = X + Y < C + S + U. 

Hence for the PD P[Z < C + S + U] of the bond issuer it holds that 


P[Z < C + S + U] < P[X < C + S + U]. 


(3.2) 







Figure 1: Illustration of Equation (3. Id): Figure 10 of Chan-Lau and Oura (2014). 


Default point Recovery upon default 2 / 


V a istfie share of encunbced assets. OC islhe overcollateralization for secioed debt, and a-A = (1+OQ-S. In 
oidefto guarantee better protection to sealed debt holders than uiseaied debt holdeis. OC {and theiefove 
a) should bet set sufficiently high so that secured debt defaults only after inseoied debt defaults (S/a < SHI 
ora>S/(SHJ))_ 

2/ Secwed debt holdeis first take min{S, a-A}. If collateral is insufficient (S>oA|, seared debt holders 
additionally receive a share of unencwnbered asset (l-a}-A in proportion to their shortfall (S-a-A). 


Secured debt 

S = a-A 


a-A + 

(l-a)-A-(S-a-A)/[U+(S- 

a-A)] 



Unsecured debt 

S+U = A 


If S< a-A, A-U 

If S>a-A, 

(l-a)-A-U/[U+( S-a-A)] 


Balance Sheet 1/ 


Asset A $100 liatoil fties $100 


collateral asset a-$A 

secured debt $S 

overcollateralization 

OC-$S 

unsecured debt $U 

unencumberd asset 
(l-a)-$A 

equity $E 


If the credit quality of the cover pool assets is much better than the credit quality of the bond issuer’s 
remaining assets, (3.2) can be a major restriction for the two-assets model because it implies that the 
issuer’s PD cannot become arbitrarily bad (see Section j.l for an illustrative example). To be more precise: 
Note that by (3.1c), the standalone (i.e. without support by the bond issuer) probability of loss of the cover 
pool is P[X < C]. Hence, if the credit quality of the cover pool is good (i.e. P[.Y < C] is small) and S+U 
is small compared to C (i.e. most of the firm’s portfolio is pledged as collateral for the covered bonds) 
then (3.2) implies that the bond issuer’s PD is also small even if the other assets are very risky. 

In any case, it is important to keep in mind that the probabilities in (3.2) may refer to different time 
horizons, depending on the purpose the model is deployed for: 

• If the model is used for calculating (say) one-year PD and LGD for the bond issuer, P [Z < C+S+U] 
refers to the issuer’s one-year PD while P[X < C + S + U] should be related to the maturity of the 
covered bonds to which the cover pool assets serve as collateral. For only after the cover pool assets 
are no longer pledged it will become clear to which extent the issuer’s unsecured creditors can take 
recourse to the cover pool assets to satisfy their claims. 

• If the model is used for calculating (say) one-year expected loss for the covered bonds, both P[Z < 
C+S + U] and P[A' < C+S + U] must refer to a one-year time horizon. For in this case a covered 
bonds loss event can only occur if the bond issuer defaults. Default of the covered bonds then is 
triggered by the value of the cover pool after one year. The covered bonds recovery value after one 
year is determined by the one-year values of the cover pool assets and the remaining assets. 

3.2. Covered bonds expected loss in the two-assets lognormal model 

Ever since the seminal paper by Merton (1974) was published, modelling of a company’s asset value 
distribution by a lognormal distribution has been very popular. To do so is convenient if not always 
realistic. Here, we follow this approach mainly because it allows fitting a model to sparse given data. 

Assume that X and Y from Section 3.1 are both lognormal variables, linked by a normal copula. Hence 

X = exp(/i + af), Y = exp(y + Trj), (3.3) 
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with fi,v £ M., a,T > 0 and (£, rf) 




for some g £ [0,1]. Since the sum of two 


lognormal random variables in general is not lognormally distributed, most of the probabilities and 
expectations of the variables listed in Section 3.1 can only be numerically evaluated. To facilitate the 
numerical calculations, in the following, we provide formulae for the evaluation of the loss probabilities 
and expected losses that do not require more numerical effort than one-dimensional integration. We need 
to distinguish the two cases g < 1 and g = 1 (when X and Y are comonotonic). 


Case g < 1. We make use of the fact that by the assumption of a normal copula linking £ and g in (3.3), 
the distribution of ?? conditional on £ = x is normal with mean gx and standard deviation yfl — g 2 . By 
applying the disintegration theorem, we then obtain the following equations for the probabilities of the 
loss events listed in Section 3.1: 


P [Z<C + S + U,Z>C + S} = E[P[r <C + S + U - X,Y >C + S - X\X]] 

l°g ( C S +£/) — /i, 

log(C+S+t/-e'‘ + ”)-(.+T g x) \ dx 


- j <p(x) $ ( los ( C+S - T e ^- (t/ + Tex) ) dx , (3.4a) 

— OO 
log (C+S)-fx 

P [Z<C + S + U,Z<C + S,X>C] = J ip(x)<f> ^°s( c +S-e^n-( , '+r QX) 

log (C) —M 


^ dx, (3.4b) 


( p{x) *1> 


-OO 

log(C + S)- 


P[Z < C + S + U, Z < C + S, X < C] 


log (C)—M 



— OO 


v{x)<P ^ G+S ~^^pY^ +Te ^ dx. 


(3.4c) 


We obtain formulae for the expected loss rates by integrating the loss variables from equations (3.1a), 
(3.1b) and (3.Id) over the respective loss events and then adding up the results. In the derivation of 
(3.5a) and (3.5b), we make use of the second equation in (2.6a). 


log (C)-^ 

CE[L c \ = J ( C - e ^ +<T3: ) y (a;) $ ^°g (c+s-e^+^-^+r ex )^ dx 

— OO 

log (C)-fj. 

_ e w +r 2 (1 — e 2 )/2 f (C - e^+ ax ) e Tex ( log (c+S-e^+ a:c )-(v+T gx) 

J C + S-e»+™ V ^ 

— OO 



dx, 


(3.5a) 
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log(C + S)-n 


SE[L s ] = (C + S)P[Z < C+S,X > C] - 


log (C)-g 


e^ +ax <p(x)<f> ( l^C+S-^'hHrgi) 


dx 


log(g)-./ 


e v+TX tp(x) { d> 


J fF) / log(C+S-e‘' +T3; )-(/J+tT g a;) 


- <E> 


log(C)-(M+°' gap 

< 7 \/i ~e 2 


dx 


log (C)-l± 


- Se l,+r2 ^-e 2 )/ 2 


V dVi-e 2 

^(z) $ fl°g(C+S-e'-+”)-(,+r g3: ) _ T dx 

V ^Vi-e 2 / 


C + S - e^+ 


+ 5P[Z<C , + 5’, X<C\, (3.5b) 


UE[L V \ = UP[Z < C + S] + (C + S + U)P[C + S < Z < C + S + U] 

log (C+S)-/x 


+ j er+^yjx) 3> ^ MC+S-e^-^+r g x) ^ 

—oo 

log(C+.S + C7)-/* 

J e M+ffa: ^(x) $ ^ los(C,+s+c/ -^ 


dx 


r~ e v+‘ Tn: )-(v+ T ex) 
- e 2 


dx 


— OO 

logCC + S)-^ 


z/+r a? 


<p(z) $ 


log(C+S— e' ,+T *) — (/J+o- gap 

V 1 -*? 2 


da; 


log(C + S + C)-i- 


j e v+T x ^x)$^c+s+u-^+-J-(»+° ( >^ dx (3 


5c) 


Case p = 1. In this case, the asset value variables X and Y are comonotonic, i.e. if one of them is 
known the value of the other is also known. This is the strongest possible type of dependence between 
two random variables. It defines a worst case scenario because high losses with the cover pool of the 
covered bonds occur exactly at the same time when also the losses with the issuer’s remaining assets are 


high. Denote by x(a) for a > 0 the unique solution of the equation 

a = e /J-+ax + e u+TX . (3.6) 

The probabilities of the loss events from Section 3.1 then can be calculated as follows 

P[Z <C+S + U,Z>C + S] = <f>(x(C + S + U )) - $(a z(C + S)), (3.7a) 

P [Z <C+S + U,Z <C + S,X>C\= max (<f>(x(C + S)) - , 0), (3.7b) 

P[Z < C + S + U, Z < C + S, X < C] = $(min (x(C + S), log( ^ ) ~ At )). (3.7c) 


Although in this case the model is basically one-dimensional, due to the non-linear structure of the loss 
variables in (3.Id), numerical integration cannot be avoided when evaluating the following equations 
(3.8a) and (3.8b) for the expected losses of the covered bonds and senior unsecured debt holders. In 
contrast, expected loss for the junior debt holders does not require much computational effort, as shown 
in (3.8c). 


CE[L c ] 


min( log( ^ ) ~ |J , x(C-TS)) 


, Ac- e» +(TX ) (C + S- e^+ CTX - e" +TX ) 

mx) --~-x- dx, 

v ’ C + S -e»+ ox 


(3.8a) 
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min( log( ^ ) ~' J , x(C+S)) 

SP[L s ] = SP[Z<C + S,X<C}-S J <p{x) c + e ^^ +(7X dx 

— oo 

+ (C + S) P[Z < C + S,X > C] - e A ‘ + °’ 2 / 2 max ($(z(C + S) - a) - $( 1 ° s( ^ ) -' x - a), 0) 

— e" +r2/2 max (<&(x(C + S) — t) - $( lo ®=^ - r),0), (3.8b) 

J7E[Lu] = t/P[A < C , + S'] + (C + S' + [/)P[Z < C + S + U,Z> C + S] 

- e^ 2 / 2 (&(x(C + S + U)~ a)- $(a;(C + S) - ct)) 

- e^ 2 / 2 ($(a;(C + 5 + U) - r) - <f>(x{C + S) - r)) . (3.8c) 

4. Calibration of the two-assets model 

The credit risk of a financial institution typically is stated in terms of the institution’s PD and LGD. 
As described in Section 2, then often a distribution of the institution’s assets value can be fitted by 
a variation of the method of moments such that the given PD and LGD are matched (see (2.1a) and 
(2.1b)). However, in order to implement a two-assets model like in Sections 3.1 and 3.2, we need the joint 
distribution of the pair (A', Y) of the values of the cover pool and of the remaining assets respectively. 

4.1. The lognormal case 

In the case of the two-assets lognormal model as specified by (3.3), we have to determine the values of 
the five parameters g, 1 a , v, r and g. It would be nice if there were a way to split the PD and LGD (or 
equivalently the PD and expected loss) associated with a financial institution that issues covered bonds 
into separate PD and LGD estimates for the cover pool and the pool of the remaining assets respectively. 
Then we could use the approach from Section 2.2 to fit lognormal marginal distributions of the cover pool 
assets and other assets respectively (i.e. the parameters g,,a, v,t) and either try to somehow estimate 
the correlation g or just choose an appropriate value for g. As there is an obvious worst case (the case 
g = 1 with comonotonic marginal distributions), choosing g might be a better approach than estimating 
it because this way the model easily can be made to err on the conservative side. 

Without further constraints, it is likely that there are many ways of splitting the issuer’s PD and LGD 
into PD and LGD parameters of the sub-portfolios, and some of these ways might even be reasonable 
from an economic perspective. But some constraints actually make sense. Given that regulatory approval 
for a covered bonds issuance is required if relief of capital requirements for the bond investors is intended, 
the cover pool is likely to be subject to a number of rules regarding its composition and riskiness. 
Therefore, it seems safe to assume that probability of loss and expected loss estimates p covel and EL cover = 
Pcover LGD coveT for the cover pool are available. 

A number of different approaches to the calibration of the cover pool asset value parameters to given 
Pcover and .ELcover are conceivable. Here we suggest calibrating to the probability of the event that the 
cover pool asset value falls below the threshold ’face value of the covered bonds times (100% plus level 
of over-collateralisation)’ and the expected value of the shortfall below the threshold. Denote by C > 0 
the face value of the covered bonds and by v > 0 the level of over-collateralisation of the covered bonds, 
as in Section 3.1. Hence we have 

Face value of the covered bonds x (100% plus level of over-collateralisation) = C (1 + u). 

For determining the parameters of the lognormal representation (3.3) of the cover pool asset value we 
make use of the observation in Remark 2.1 about the relationship between PD, threshold and recovery 
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rate on the one hand and PD, quantile and ES on the other hand. Once we have expressed quantile and 
ES of the asset value in terms of the threshold C (1 + v) and EL cover = p cover (1 — RR C over), we can use 
results from Section 2.2 to infer the values of /i and a. 

Let A' denote the asset value of the cover pool. As suggested in Remark 2.1, we then have C (1 + v) = 

tfpcoverPO and 


ES V (X) 

I'cover \ / 

C (1 + v) 


RRc 


ES Pcover { X) = C{ 1+v) 


Pcover EL C 
P cover 


(2.6a) and (2.6b) now imply the following equation system for p and a: 


/i = log((l + v) C) - (7 $ 1 (pcover )? _ 

0 = (.Pcover) ( Pcover EL cover) (Pcover) & /2). 

The factor l+v in the first row of (4.1) reflects the fact that the risk characteristics are provided for the 
whole cover pool, not only for that part of it which has exactly the value of the covered bonds. Recall that 
according to Proposition 2.3, (4.1) has a unique solution (p,cr) eRx (0, oo). Parameters p and er of the 
cover pool asset value hence can be assumed to be known. This is not the case for the parameters v and 
r of the distribution of the value of the remaining assets in the bond issuer’s portfolio because typically 
no detailed information is publicly available about the structure and riskiness of a bank’s portfolio unless 
the assets are pledged as collateral for borrowing. 

Nonetheless, we can try and determine implied values for v and r as follows. Denote by Pi SSU er the 
issuer’s PD, by LGD issuer the issuer’s LGD 6 and by Assuer the issuer’s total amount of outstanding debt 
(including the covered bonds issued). In the setting of Section 3.1, it would hold that Assuer = C + S + U. 
As in Section 3.2 we then represent pi SS uer and LGD lssuel as functions of the parameters p, a , u, t and g. 

In case of g < 1, with A' and Y as in (3.3), we obtain as the first equation 


Pissuer — P[A + Y < -A ssuer 


lo S( D isauer)-^ 



—OO 


(4.2a) 


For the second equation, we obtain 


-A SSU er Pissuer (1 EGD[ ssuer 'j — E [(A + L ) 1 { X-r>'< D; ssu( . r } ] 

log(Dissuer)-M 


J C fJ,+aX (f(x) T ^ 1 * * °g( D i^'r°r~^ + ' TX )~l iy + ' r g X ) ^ 


dx 


— OO 


e v + r2 (!- e 2 )/ 2 
In case of g = 1, the first equation is 


j e Tex v(x) $ ^ lo s(^suer-e^+-D-(.+r ex) _ T dx , ( 4 . 2 b) 


Pissuer 


^ (-^issuer ) ) •> 


(4.3a) 


6 For the sake of simplicity, we assume that we are given an LGD for the issuer’s entire debt. In practice, it might be more 

likely to have an LGD for the issuer’s senior unsecured debt. This case could be dealt with by replacing the right-hand 

side of (4.2b) with the right-hand side of (3.5b) in case of g < 1 and the right-hand side of (4.3b) with the right-hand 

side of (3.8b) in case of g = 1. 
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where for a > 0 the function x(a) is defined as the solution of (3.6). The second equation is given by 

-^issuer Pissuer (1 LGD[ ssneY ^ — C^ ^ ^ (*^(-^issuer) a) +e v+T 2 / 2 <S>(x(D issuer) r). ( 4 . 3 b) 

Assuming that parameters p, er and g are known, we have to solve the equation system (4.2a) and (4.2b) 
(or (4.3a) and (4.3b) respectively) for v and r in order to fully parameterise the lognormal two-assets 
model from Section 3.2. Unfortunately, there is no obvious way to transform the equation system (4.2a) 
and (4.2b) such that variables v and r are separated. This is different for equation system (4.3a) and 
(4.3b) in the comonotonic case g = 1. We can therefore state the following result on the solutions of the 
system. The proposition shows that there are combinations of the input data p cover , EL covel , G, v, pi ssuer , 
-bG-D issuer; and Assuer suc h that the equation system (4.3a) and (4.3b) for v and r has no solution at all. 


Proposition 4.1 Assume that Assuer > 0, 0 < pi ssuer < 1, 0 < LGAssuer < 1, p € R. and a > 0 are 
known and kept fix. Then there exists a solution ( v , r) £ R x (0, oo) of the equation system (4.3a) and 
(4.3b) if and only if it holds that 


0 < Pis 


< 


f log( Assuer ) P 


Pissuer e A ‘ +<T $ "(Pissuer) _ e M+<U/2 l(p. ssuer ) _ cr) 


Pissuer h)\. 


< LG Du 


< 1 - 


and (4.4a) 

6 m+- 2 /2$($-i (p issuer )-a) 


Pissuer -A; 


(4.4b) 


If there is a solution {v, r) elx (0,oo) of (4.3a) and (4.3b) it is unique. 


Proof. For the sake of a more concise notation, in this proof we drop the index ’issuer’. Observe that for 
any finite (v,t) the solution x(D) of the equation 


must satisfy 


£) — e ^+<T x(-D) _|_ e v+rx(D) 

X (D) < ! 2 S<£Uii. 


(4.5) 


Comonotonicity (g = 1) implies that the random variables X and Y from (3.3) can be represented as 
X = e' i+ °A, Y = e tJ+T £ for some standard normal variable £. Hence it follows that 


p = P[e'* +<7 * + e u+Ti <D}= P[f < x{D)} < P 


£< 


log (D) - g 


= $ 


^ log {D)-p 


This shows that (4.4a) is necessary for the existence of a solution. 

Assume now that p satisfies (4.4a) and that (v,t) is a solution of (4.3a) and (4.3b). Then (4.5) can be 
solved for v: 

v = log (D-e fJ ' +ax{D) ) ~tx(D). 

Substitute this for v in the term e ! ' +r ’"/ 2 $(x(D) — r) of (4.3b) to obtain the term 

g(r) d = (£> — e / i +°' x C- D )) e r2 / 2 ~ TX ( D ) $(a"(D) — t). 

The derivative of g with respect to r is 

g'(r) = (.D — e ll+ax( ' D ^) e T l 2 ~ TX i D ) |$(a;(T)) — r) (r — x(D)) — <p(r — x(D))J. 
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Figure 2: Illustration of Proposition 4.1. Ranges of bond issuer’s PD and LGD for two set of parameters 
for the cover pool. The smaller range (doubly shaded) is related to cover pool PD and LGD of 
0.05% and 30% respectively. The larger range (singly shaded) is related to cover pool PD and 
LGD of 0.5% and 50% respectively. See Example 4.2 for further explanation. 



By (2.3), if follows that g'(r) < 0 for all r > 0 such that g is a strictly descreasing function. Hence g is 
an ’onto’ mapping from (0, oo) to (lim^oo g(r), lini T _>o c/(t)) ■ It is easy to see (using again (2.3)) that 

lim T ^ o5 (r) = (D-e^^)p, 
lim-r-^oo g{r) = 0. 

By rearranging (4.3b) such that LGD is the only term on one side of the equation, it follows that (4.4b) 
is necessary for the existence of a solution of (4.3a) and (4.3b). Sufficiency follows from the fact that g is 
’onto’. Uniqueness of the solution is a consequence of the fact that g is strictly decreasing. □ 

Example 4.2 To illustrate the consequences of Proposition f.l, in Figure 2 we plot the possible PD-LGD 
ranges of the bond issuer for two sets of parameters for the cover pool (see Section 3.1 for information 
on the meaning of the parameters): 

• C = 0.3, S = 0.6, U = 0.1, v = 0.2, 

. {pcover , LGD cover ) € {(0.0005,0.3), (0.005,0.5)}, 

• g= 1 (comonotonic case). 

Reading example for Figure 2: PD = 2% for the bond issuer cannot be realised with either of the two 
parameter sets for the cover pool. PD = 0.8%, LGD = 40% for the bond issuer can be realised with 
(Pcover, LGD covei ) = (0.005,0.5) but not with (p cover , LGD cover ) = (0.0005,0.3). 

Intuitively, the restricted PD-LGD range for the bond issuer from Proposition 4.1 and Example 4.2 is 
caused by the following phenomenon: The input parameters p cover , EL cover , C , and v for the distribution 
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of the cover pool asset value X may force X to have a large variance while at the same time the input 
parameters pi ssuer , LGDi ssuer , and Z?i ssuer for the distribution of the issuer’s asset value Z = X + Y may 
force Z to have a small variance. As a consequence, it may happen that model (3.3) is not able to fit the 
given risk characteristics. But recall the comments from Remark 3.1 regarding the appropriate choice of 
time horizons for the risk characteristics. To some extent, with the right time horizons, the problem of 
not having an exactly fitting two-assets model thus might be mitigated. 

Due to the complex structure of the equation system (4.2a) and (4.2b) it would be hard to try and analyse 
the particulars of the question when exactly there is no solution to the system. To better illustrate the 
cause of the problem, in the rest of the section we instead analyse a toy model based on two jointly 
normal distributed asset value variables. 


4.2. Model existence in the case of two jointly normal asset values 


In this example, we leave all assumptions from Sections 3.2 and 4 unchanged but replace (3.3) with 

X = fi + a£, Y = v + tt], (4.6) 

"o\ A A' 


where we have g, v £ R, a, r > 0 and (£, 77 ) ~ Af 


for some g £ [0,1]. This is not a 


y 0J 1 \g 1 , 

satisfactory assumption for modelling asset values because X and Y may take negative values. But it 
is helpful for the purpose of illustrating potential calibration problems because the sum of two jointly 
normal random variables is again normal. As a consequence, in this case, in contrast to the lognormal case 
of Equations (4.2a) and (4.2b) (or (4.3a) and (4.3b)) we obtain a closed-form solution for the calibration 
equations. 

Under (4.6) the equations for the parameters g and a are given by (cf. Section 2.1) 


P = (1 +v)C-a$ 1 (p cover ); 

(1 4- v) C EL cove] 


a = 


(4.7) 


Pcover d* (p cover ) + <p(d> 1 (p C over)) 

The equation system (4.2a) and (4.2b) for v and r is replaced by 

^ = Assuer p U T (Pissuer): 

t = -vA 2 - (1 - g 2 )cr 2 - a g, 


iP = 


-^issuer Pissuer bCi Dj ss)]f T 


(4.8) 


Pissuer d* (pissuer) A pA* (Pissuer)) 

With regard to (4.8), it is easy to see that 

r > 0 ip > a. 


(4.9) 


The parameter ip is actually the implied standard deviation of the bond issuer’s total asset value. As 
observed above, condition (4.9) is natural for the existence of a two-assets model with PD and LGD 
constraints for the cover pool and the issuer’s entire portfolio. However, while in the case of assumption 
(4.6) the condition is indeed necessary and sufficient for the existence of a model based on a joint normal 
distribution, it is unclear whether a similar statement holds true in the case of assumption (3.3). Note 
that even if (4.9) is satisfied the resulting model still might not make economic sense because it could 
occur that v < 0. 
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5. The lognormal one-asset model 


As shown in Section 4, it may be impossible to exactly fit a two-assets model for the covered bonds 
calculations. Therefore, in this section we consider how to implement a lognormal one-asset model and 
explore how to make it sensitive with regard to differences in the riskiness of cover pool and remaining 
assets. 

Lognormal one-asset model. The uncertain future value A of the total assets of the covered bonds 
issuer (including the cover pool assets) is given by 

A = exp (re + V’Ot (5-1) 

where £ is a standard normal random variable. The value of the cover pool assets is e A for a constant 
£ £ (0,1], The number e is called asset encumbrance ratio. 

In this section, basically we again assume the setting of Section 3.1, with X = e A and Y = (1 — e) A: 

• The calibration of the parameters k and ip of the lognormal one-asset model then is straightforward 
(see the comments on the derivation of (4.1)) if the issuer’s PD pi ssuer and LGD LGD\ ssue r are 
known: 

« = log(C + S + U) - 'Ip $ -1 (pissuer), ^ 

0 = (.Pissuer) 'P') (Pissuer -S'-^issuer) (Piusser) /2) , 

with EL[ ssuer = Pi SS uer AGD; ssuer . By Proposition 2.3 the equation system (5.2) has a unique solution 

(K,ip). 

• In principle, the encumbrance ratio e could be determined as 

^ = ( 1 + t; ) AT- ( 5 - 3 ) 

^*-0 

where A 0 > C + S + U are the issuer’s total assets at time 0 (today), v is the level of over¬ 
collateralisation of the covered bonds and C is the face value of the covered bonds. 

In practice, the value of Aq might not be accurately known, for a number of possible reasons. In particular, 
under an amortised cost accounting regime the sum of the asset values might be considered a too unreliable 
and potentially too conservative estimate of their true total value. That is why, in this section, we explore 
how to infer e implicitly from the level of over-collateralisation v and the expected loss of the cover pool. 

5.1. Adjusting the asset encumbrance ratio 

Assume that in addition to the debt amounts C, S , and U and level of over-collateralisation v as in 
Section 3.1 and parameters n, ip (calibrated according to (5.2)) a value EL covel of expected loss for the 
cover pool is given. 

The idea is to find e such that 

(1 + v)CEL cover = E[((l + v) C — e A) l{ e a<(i+v) c}] ■ (5-4) 

Actually, (5.4) is equivalent to the second equation of (4.1). Recall that we represent the value of the 
portfolio as A = e K+1 ^ with a standard normal random variable £. With •& = K + log(e) therefore we may 
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rewrite (5.4) as 


(1 + v)CEL covei = E[((l + v)C-e^t)l {e <, + ^ <{1+v)c} ] 

= (1+ V )C* (^ l0g((1 + J } C) ~ ^ - e^ 2 / 2 $ ^ 0g((1 + J } C) ~ * 

d = m. 



Based on (5.5), it can be shown that (5.4) always has a unique solution e: 


Proposition 5.1 There is exactly one id G K such that /(d) = (1 + v) C EL cover , with /(d) defined hy 
(5.5). 


Proof. Obviously, we have lim^-oo /(d) = (1 + v) C. The limit behaviour lim#-**, /(d) = 0 is less 
obvious but follows from (2.3). Since 0 < EL covel < 1, therefore by continuity of / there is a solution d 
to /(d) = (1 + v) CEL covei . Since di < d 2 implies < (1 + v) C} C {e ,Jl+ ^£ < (1 + v) C } we 

obtain for di < d 2 

E[((l +v)C- e*+^)l {e * 1+ , t<(1+t , )c} ] > E[((l + «) (7 - e*' + *e)l {e * a+ + t<1+v)C} ] 

> E[((l + v) C - e^+^S)l {e , 2+ « <(1+ , )c} ]. 

This implies uniqueness of the solution of /(d) = (1 + v) C EL covei . □ 

By Proposition 5.1, we can define a lognormal one-asset model for covered bonds expected loss in the 
following way: 


Definition 5.2 Assume that the issuer's PD pi SSU er, the issuer’s expected loss ELi ssuel , the cover pool 
expected loss EL cover , the debt face values C for the covered bonds, S for the senior unsecured debt, U 
for the junior unsecured debt, and the level of over-collateralisation v of the covered bonds are given. The 
adjusted lognormal one-asset model then is specified by the following two calibration steps: 

• Parametrise the asset value A = e K+ ^ with f ~ Af( 0 , 1), k £ ffi. and ip > 0 by solving (for k and 
ip) the equation system (5.2). 

• Determine the adjusted encumbrance ratio e by solving (5.5) for the unique solution d and setting 
e = min(l,e’ ? - K ). 

The 100% cap of e in the second calibration step of Definition 5.2 is needed because while e is always 
positive there is no guarantee that it is smaller than 1, in particular if the values of EL cover and EL, ssuer 
are very different. A value of e greater than 100% would not make economic sense because then the value 
of the non-cover-pool assets in the issuer’s portfolio would be negative. Actually, in a way similar to 
Proposition 4.1, we can state which values can be exactly realised for EL covel for given characteristics of 
the issuer’s portfolio . The following observation is a consequence of the fact that the function /(d) in 
(5.5) is a decreasing function of the parameter d (as shown in the proof of Proposition 5.1). 


Remark 5.3 In the setting of Definition 5.2, the minimum expected loss of the cover pool that can be 
matched with the adjusted lognormal one-asset model is 


EL 


cover, min 


/log((1 + v) C) - ti\ 

\ * ) 


gK+i/d /2 /log((l + v) C) — K 
(1 + v) C ^ y Ip 



How different is the adjusted lognormal one-asset model from the lognormal two-assets model as discussed 
in Sections 3.2 and 4.1? The following proposition is useful for figuring this out. 
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Proposition 5.4 Consider the following two bivariate random vectors: 

• (,e A, (1 — e) A) with 0 < e < 1, A = e K+ ^ < ’, ( ~ 7V(0,1), k £ M, ip > 0. Call this Model 1. 

. (X,Y) withX = e»+°t, Y = e v +^, Q ~^((q)»Q j))» a,r >0, ? € [-1,1], 

Call this Model 2. 

Model 1 and Model 2 are equivalent in the sense that the respective random vectors have the same distri¬ 
butions if and only if 


ip = t = a, g = 1 , e = 


ev + e' 


r. « = log(e M + e"). 


(5.6) 


Proof. Assume first that (5.6) is satisfied. By g = 1 and a = r we then have almost surely that 

(X,Y) = (e" 


— e v + cr ^\ 


Making use of ip = a and of the representations of e and k in (5.6) we obtain 


= £ e K+ ^, 




= (l- £ ) 


0 «+VA 


Since £ ~ A/"(0,1), this implies the ’if’ part of Proposition 5.4. 

Assume now that the bivariate distributions implied by Model 1 and Model 2 respectively are equal. 
Then in particular, the distributions of log(e A) = log(e) + k + ip£ and log(X) = p + a tp are equal. 
Since both £ and ^ are standard normal, it follows that /i = log(e) + n and cr = ip. The same rationale 
applied to (1 — e) A and Y gives v = log(l — e) + k and t = ip. Hence we have ip = t = a and 
g. — log(e) = v — log(l — e) which gives e = el f +e „ ■ This implies n = // —log(e) = log(e M +e y ). By ip > 0 the 
components of (e A, (1 — e) A) are comonotonic. Hence also X and Y are comonotonic which implies that 
log(X) = g + a f and log(H) = v + r ?? are comonotonic. However, two jointly normal random variables 
are only comonotonic if their correlation is 1. This completes the proof of ’only if’. □ 

By Proposition 5.4, for homogeneous portfolios, the adjusted lognormal one-asset model of Definition 5.2 
is a comonotonic special case of the lognormal two-assets model from Section 3.2. To formally show this, 
we consider the following type of lognormal two-assets model. 

Definition 5.5 A lognormal two-assets model as specified by (3.3) is called margins-calibrated if it is 
calibrated in the following way: 

• Debt face values C for covered bonds, S for senior unsecured debt and U for junior debt are given. 

• A level v of over-collateralisation is given. 

• A correlation parameter g is chosen on the basis of expert judgment or a separate calibration exercise. 

• Parameters fi and a are determined by solving the equation system (4.1). 

• Parameters v and t are determined by solving the equation system 

V = log(5 + U - V C) - T $ _1 (Pother), 


0 — (Pother) y) (Pother Another) CXp (t 4* (Pother) T /2) • 


(5.7) 


where p 0 thei and EL ot h er denote the stand-alone probability of loss and the expected loss respectively 
of the issuer's portfolio without the cover pool. 

The term —vC is a consequence of having the term 1 + v in (4.1) to reflect the over-collateralisation 
of the cover pool. The term is a correction needed to make sure that the total debt is C + S + U. The 
calibration approach in Definition 5.5 is computationally less demanding than the calibration approach 
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from Section 4.1 and always feasible. But as mentioned in Section 4.1, often no estimates of Pother and 
E Mother might be available as input to (5.7). Nonetheless, for comparison with the adjusted lognormal one- 
asset model and illustrative purposes, the margins-calibrated two-assets model is helpful as demonstrated 
in the following remark and Section 6 respectively. 


Remark 5.6 In the comonotonic case (g = 1) of Definition 5.5, if in addition Pother = Pcovei and 
EL 0 ther = EL covel (i.e. the issuer’s portfolio is homogeneous in terms of risk) then by Proposition 2.3 
if follows that cr = r. Hence by Proposition 5-4 the two-assets model can be represented as a lognormal 
one-asset model. This one-asset model actually is an adjusted lognormal one-asset model in the sense of 
Definition 5.2 with pi ssuer = Pcover and EL[ ssuel = EL covel . This follows from the observations that 

• with e = the calibration equation systems (5.2) for the one-asset model and (4.1) and (5.7) 

for the two-assets model are all equivalent, 


• by (5.6) we have = ee K , and 

• (4.1) implies (5.4) with e A replaced by e A1+<T ^ where £ is standard normal. 


5.2. Expected loss for the adjusted lognormal one-asset model 


In the case of a lognormal one-asset model like (5.1), with encumbrance ratio e calibrated according to 
(5.3) or as in Definition 5.2, it is straightforward to derive the expected loss formulae for the loss variables 
defined by (3.1a), (3.1b) and (3.Id) (with X = e A and Y = (1 — e) A). For the readers’ convenience, we 
list these formulae in the following. 

• Probabilities of the loss events: 


P [A<C + S + U,A>C + S\ = Q 
P[A<C + S,eA>C] = $ 
P[A<C + S, eA<C] = $ 


- $ 


^ log((7 + 5 + l7))-/x 

/ log(max(C 1 /£, C + S)) - p \ _ ^ / log [C/e) - g 


/ log (C + S)-p 


log(min(C/e, C + S)) — g' 


(5.8a) 

(5.8b) 

(5.8c) 


• Expected loss for the different types of debt (covered bonds, senior unsecured and junior unsecured): 


CE[L c \ 


SE[L s ] 


UE[L V ] 


log(min (C/e,C + S))-ia 


, (C — se ll+(TX ) (C + S— e M+<Ta: ) 

n x ) -e;—- dx , 


C + S~ee»+° 


(5.9a) 


—oo 

log(min(C/£,C + S))-^ 


S 


c + S- ee»+™ dx + (C + S)P[A<C + S,eA>C] 


_ e fj,+a 2 /2 | $ ^ log(max(C/g,C , + S)) - \i _ ^ $ ^ log((7/e) - g _ ^ j 


(5.9b) 


UP[A < C + S] + (C + S + U) P[A < C + S + U, A > C + S] 

- ( W + s + ?) - p _ „ ( i.g (c + s) - p _ . (5 , 9c) 
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Table 1: Illustration of expected loss (as % of exposure) for covered bonds, senior unsecured debt and junior debt 
as function of the correlation (as %) between the cover pool logarithmic asset value and the logarithm 
of the value of the remaining debt. Approach and parameters are explained in Section 6.1. 


Correlation g 

0.0 

30.0 

60.0 

90.0 

100.0 

Covered bonds EL 

0.002 

0.014 

0.057 

0.173 

0.257 

Senior unsecured EL 

0.007 

0.039 

0.132 

0.345 

0.465 

Junior EL 

0.020 

0.097 

0.294 

0.711 

0.943 

All EL 

0.007 

0.038 

0.126 

0.330 

0.450 


6. Numerical examples 

Due to the total lack of credit loss data related to covered bonds, we cannot directly test the models 
described in Sections 3.2, 4.1 and 5.1. What we can do is checking model results for economic plausibility 
and comparing the margins-calibrated two-assets and adjusted one-asset models. Hence, in this section, 
we present some results to illustrate how the models work in the following three scenarios: 

• Increasing correlation between the values of the cover pool assets and the issuer’s other assets 
respectively. 

• Increasing asset encumbrance by increasing the proportion of the total debt consisting of covered 
bonds. 

• Using the adjusted lognormal one-asset model instead of the lognormal two-assets model. 

6.1. Impact of dependence between cover pool and remaining assets 

We use the margins-calibrated lognormal two-assets model as described in Definition 5.5 to illustrate the 
impact of correlation on the expected loss of covered bonds (senior secured debt), senior unsecured debt 
and junior (unsecured) debt respectively. The input parameters are as follows: 

• C = 0.3, S = 0.6, U = 0.1, v = 0.2. 

• Pcover ~ Pother = 0.01, LGD cover — LGD 0 th er — 0.45. 

• g £ {0, 0.3,0.6, 0.9,1.0}. This generates a spectrum of dependence increasing from independence 
(g = 0) to comonotonicity (g = 1). 

The CGFS recently found in a survey of 60 European banks 28.5% median asset encumbrance in the 
banks’ balance sheets (CGFS, 2013, Graph 3). This finding was the motivation for the choice of the debt 
amounts in the debt hierarchy of this example, together with the consideration that the proportion of 
junior debt in general is small compared to the proportion of senior debt. 

The results of the calculation are shown in Table 1. The impact of increasing correlation on the expected 
loss is very strong across all debt types but strongest for the covered bonds (senior secured debt). This 
behaviour is intuitive because strong dependence implies a larger likelihood of a big part of the issuer’s 
assets being erased. In the following two examples, we consider the comonotonic case (g = 1) only, thus 
adopting a worst case perspective that, in particular, might be appropriate for stress testing purposes. 

6.2. Impact of asset encumbrance 

Here we use the adjusted lognormal one-asset model as described in Definition 5.2 to illustrate the impact 
of asset encumbrance on the expected loss of the senior unsecured debt and the junior (unsecured) debt 
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Table 2: Illustration of expected loss (as % of exposure) for covered bonds, senior unsecured debt and junior debt 
as function of the covered bonds exposure (as % of total debt). Approach and parameters are explained 
in Section 6.2. 


Covered bonds exposure C 

0.0 

10.0 

20.0 

30.0 

40.0 

50.0 

60.0 

70.0 

80.0 

Adjusted encumbrance ratio 

0.0 

12.0 

24.0 

36.0 

48.0 

60.0 

72.0 

84.0 

96.0 

Covered bonds EL 

NA 

0.237 

0.246 

0.257 

0.269 

0.283 

0.300 

0.320 

0.347 

Senior unsecured EL 

0.395 

0.415 

0.438 

0.465 

0.496 

0.536 

0.587 

0.658 

0.777 

Junior EL 

0.943 

0.943 

0.943 

0.943 

0.943 

0.943 

0.943 

0.943 

0.943 


respectively. The input parameters are as follows: 

• Pissuer = 0.01, LGD[ ssuer = 0.45. 

• EL covei = 0.0045. 

• C G {0,0.1, 0.2,..., 0.8}, S = 0.9 — C, U = 0.1, v = 0.2. This generates a spectrum of asset 
encumbrance increasing from none (C = 0) to strong (C = 0.8). 

Note that according to Remark 5.6, this specific adjusted one-asset model (with ELi asuel = EL cover ) is 
equivalent to the margins-calibrated model with p covel = Pother and LGD cover = LGD ot ^ el . Moreover, 
the adjusted asset encumbrance ratio in this case is e = 

The results of the calculation are shown in Table 2. Due to the subordination of the junior debt when 
it comes to distributing the recoveries, the expected loss for the junior debt does not depend on the 
proportions of secured and unsecured senior debt. There is, however, a significant dependence of the 
senior unsecured expected loss on the asset encumbrance ratio. The dependence becomes the stronger 
the more the asset encumbrance ratio approaches 100%. But also the per unit expected loss of the covered 
bonds increases with increasing asset encumbrance ratio. This is caused by the fact that the buffer for 
the covered bonds provided by lower ranking and unsecured debt shrinks with the increasing proportion 
of the total debt consisting of the covered bonds. Note that nonetheless in Table 2 irrespective of the 
changes in the covered bonds exposure the expected loss of the issuer’s entire portfolio remains constant 
at 0.45%. 


6.3. Lognormal two-assets vs adjusted lognormal one-asset 

Remark 5.6 presents sufficient conditions to coincide for the margins-calibrated lognormal two-assets 
model as in Definition 5.5 and the adjusted lognormal one-asset model as in Definition 5.2. In this section, 
we show that there is no longer equivalence if the assumption Pother = Pcover and £X ot her = EL covel is 
violated, i.e. if the cover pool and the pool of the issuer’s remaining assets are no longer homogeneous in 
terms of risk. 

With regard to the debt exposures and the level of over-collateralisation, we use the same input parameters 
as in Section 6.1, i.e. 

C = 0.3, S'= 0.6, [7 = 0.1, v = 0.2. 

For the two-assets model according to Definition 5.5, we use use the following parameters: 

• g = 1, i.e. we consider the comonotonic case. 

• Pother “ 0.01, TG'D 0 ther = 0.45. 

• We keep EL covel fix at 0.003, but choose LGD covel £ {0.3,0.45,0.6} and determine p c OV er by 

FIT, 

Y) — ■ L/ - L 'cover 

Fcover t,C1 D 

Js cover 
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Table 3: Results of the margins-calibrated lognormal two-assets model compared with results of the adjusted 
lognormal one-asset model. For model calibration and input parameters see Section 6.3. All numbers 
are percentages. 


LGD cover 

30 

45 

60 

V cover 

1.000 

0.667 

0.500 

Pissuer 

1.000 

0.848 

0.655 

ETr 

- L/ -^issuer 

0.396 

0.382 

0.348 


2 assets 

1 asset 

2 assets 

1 asset 

2 assets 

1 asset 

Covered bonds EL 

0.147 

0.168 

0.189 

0.188 

0.216 

0.211 

Senior unsecured EL 

0.431 

0.421 

0.409 

0.409 

0.368 

0.371 

Junior EL 

0.931 

0.929 

0.799 

0.799 

0.629 

0.628 


The corresponding one-asset model according to Definition 5.2 is calibrated as follows: 

• We determine the parameters fi, a, v and r for the two-assets model as described in Definition 5.5. 

• We then compute pi ssuer as the sum of the probabilities given by (3.7a), (3.7b) and (3.7c). Similarly, 
we compute ELi ssuel as the sum of the right-hand sides of (3.8a), (3.8b) and (3.8c), divided by 

c + s + u. 

• As input parameters for the calibration of the one-asset model as described in Definition 5.2 we use 
^issuer and ELi ssuel to determine k and ip, and EL covel and v to determine e. 

Table 3 shows that as long the risk characteristics of the issuer’s entire portfolio are identical (in terms 
of PD and expected loss) as input to the two models and the difference in expected loss of the cover pool 
and the pool of the remaining assets is not too big, the resulting expected loss values for the different 
debt types according to the two-assets model and the one-asset model are quite close. Material differences 
in the results by the two models may be incurred especially if the cap of the asset encumbrance ratio as 
suggested in Definition 5.2 kicks in. This may happen if the expected loss of the issuer’s entire portfolio 
is much larger than the expected loss of the cover pool. Then it is possible that the adjusted one-asset 
model outputs an expected loss value for the covered bonds which is larger than the expected loss of the 
cover pool - a quite unintuitive result. 

Another interesting observation from Table 3 is that the shape of the cover pool asset value distribution (as 
determined by its probability of loss and loss-given-default rate) has some impact on both the portfolio¬ 
wide expected loss and the expected losses of the debt types. In particular, the more skewed the cover 
pool asset value distribution becomes, the larger the covered bonds expected loss grows. 


7. Conclusions 

Due to the growing importance of covered bonds as a funding tool for banks and as investment for 
risk-averse investors, quantitative approaches to expected loss for covered bonds and the impact of as¬ 
set encumbrance by covered bonds on expected loss for senior unsecured debt are of high interest for 
practitioners and researchers. 

In this paper, we have studied a natural lognormal two-assets model for the expected loss of covered 
bonds. It has turned out, however, that in some circumstances exact calibration of this model is impos¬ 
sible. Therefore, we have suggested a simple and computationally efficient adjusted lognormal one-asset 
approach to these issues. Input data requirements by this approach are light: PD and LGD of the covered 
bonds issuer, expected loss of the cover pool as well as the distribution of debt types in the issuer’s 
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portfolio and the level of over-collateralisation of the covered bonds. The model therefore can be applied 
in situations like stress testing where little information is available and many bonds have to be evaluated 
in short time. 

We have demonstrated that for homogeneous portfolios the one-asset approach is a special case of the 
more natural two-assets approach whose calibration, however, is difficult and sometimes not feasible at 
all. Backtesting of the two models discussed is impossible because never any covered bonds defaults 
were observed. Nonetheless, we have presented numerical examples that suggest that the results by the 
proposed ’adjusted lognormal one-asset model’ concur with intuition as long as the risk levels of the 
issuer’s entire portfolio and of the cover pool of the covered bonds are not materially different. 
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A. Appendix: Mean-variance matching 


Fitting a two-parameter distribution to given values of mean and variance is a common approach in 
finance, in particular when little is known about the target distribution otherwise. This approach is 
referred to as mean-variance matching and is a special case of moment matching or the method of 
moments. 

We start the description of mean-variance matching by revisiting the location-scale family associated 
with some random variable X in the sense of Definition 2.2. Obviously, if X has a positive and finite 
variance, for each pre-defined pair (/x, a) £ R x (0, oo) there is exactly one distribution in the location-scale 
distribution family of X with mean /i and standard deviation a. To identify the parameters m and s of 
this distribution, we have to solve the equation system 


Hence we have 


fj. = E[m + sA'] = m + sE[A] 
er 2 = var[m + sA] = s 2 var[A]. 


s = 


\/var[A]' 


171 = 11 


\/var[A] 


E[A], 


(A.la) 


(A.lb) 


Quite often in (A.lb), the variable X will be standardised such that it follows that m = /x and s = a. 
In particular, this applies to the case where X is standard normal. Occasionally, it may be inconvenient 
to standardise X as in case of the t-distribution family with fixed degree of freedom n > 2. Then the 
’standard’ distribution of the family would be characterised by E[X] = 0 and var[A] = " 2 > 1. 


Exponential transformation of location-scale families. As discussed in Section 2.1, for generating two- 
parameter families of distributions on the positive real half-axis is convenient to take the exponential of 
a location-scale family in the sense of Definition 2.2. In general, in this case there is no longer a closed- 
form solution to the mean-variance fitting problem. In the important special case where A is a standard 
normal variable, (2.5) just gives the family of lognormal distributions. It is common knowledge that the 
mean-variance fitting problem 

/x = E[exp(m + s A)] = exp(m + s 2 /2), 
a 2 = var[exp(m + s A)] = (exp(s 2 ) — l) exp(2 m + s 2 ), 


with n > 0 and a > 0 then has the following solution (see, e.g., Wikipedia contributors, 2015b): 


s = 


m = 



(A.2b) 
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Normal transformation of location-scale families. Two-parameter families of distributions on the unit 
interval can be generated by applying a mapping of R onto (0,1) to a location-scale distribution family 
in the sense of Definition 2.2. In principle, each distribution function of a distribution on the entire 
real axis (or link function) would do the job. However, for arbitrary combinations of link function and 
location-scale distribution families, in general there are no closed-form solutions for the mean and the 
variance. 

The case of the normal location-scale family with standard normal link function (denoted by $ in the 
following) is a notable exception. We are then studying distributions of the form 

Y = $(m + sX), (A.3) 


where X denotes a standard normal random variable and we assume m £ R and s > 0. With the 
alternative parametrisation 

$ ^(p) I ~Q 

m = - , s = J- -, 0<p<l,0<p<l, 

V 1 ~ Q V 1 “ e 

the distribution of Y given by (A.3) is called Vasicek distribution by some authors. It is easy to show 
(see, e.g., Meyer, 2009) that 

E[f(m + sX)] = $(m), 

2 (A.4) 

var[<b(?n + s A')] = <b 2 (m, m; ) - 4>(m) 2 , 

with $ 2 (-, -;r) denoting the bivariate standard normal distribution with correlation r. Note that for any 
random variable Z with values in [0,1] and P[0 < Z < 1] > 0 it holds that 


0 < E [Z] <1 and 0 < var[Z] < E [Z] (1 - E [Z]). (A.5) 

As a consequence of (A.4) and (A.5), it can be shown that for any pair (/x, cr) with 0 < /x < 1 and 
0<cr 2 </x(l — /x) there is exactly one solution (m, s ) of the equation system 


H = $(m), 

cr 2 = $ 2 (m,m; - 4>(m) 2 . 


(A.6) 


This solution is given by m = $ 2 (/x) and s = 
equation 



where r is the unique solution in (0,1) of the 


0 = ( C > 2 (to, to; r) — <b(m) 2 — cr 2 . 


Hence there is no fully explicit solution of the mean-variance fit problem for the Vasicek distribution. 
That is why despite the economic credentials of the Vasicek distribution (Meyer, 2009) in practice often 
the two-parameter family of the Beta distributions is used for mean-variance fit in the unit interval. 

The Beta distribution with parameters a > 0 and /J > 0 is specified by its density 


K,p{ x ) 


r(q+ff) 

i»r(/3) 

0, 


x a_1 (1 


— x)P 1 dt , 


0 < x < 1, 
x < 0 or x > 1, 


where E denotes the Gamma-function T(a) = / 0 °° t a_1 e _t dt for a > 0. If /x and cr with 0 < fi < 1 and 
0 < a 2 < /x (1 — /x) are given and Y denotes a Beta-distributed random variable, the parameters a and /3 
of the distribution of Y can be determined by solving the equation system 


M = E[V] 


a 

a + /3 ’ 


cr 2 = var[V] 


a/3 

(a + ft) 2 (a + P + 1) ’ 


(A.7a) 
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for a and /3. It turns out (see, e.g., Wikipedia contributors, 2015a) that there is a closed-from solution 
given by 


a 

P 




m(i - m) 


- 1 


(i-M) 


m(i - M) 



(A.7b) 
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